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Abstract 

We analyze the dynamical Casimir effect for a massless scalar field confined 
between two concentric spherical shells which impose on the field mixed bound- 
ary conditions. We thus complement a previous result [Phys. Rev. A 78, 032521 
(2008)], where the same problem was considered but in that case the field was sub- 
mitted to a Dirichlet boundary condition in both moving spherical shells. A general 
expression for the average number of created particles is deduced for an arbitrary 
law of radial motion of the spherical shells. This expression is then applied to har- 
monic oscillations of the shells and the number of created particles is analyzed and 
compared with the results obtained under Dirichlet-Dirichlet boundary conditions. 
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1 Introduction 



Since the pioneering paper published by Moore in 1970 |[T] and the contributions by 
Fulling and Davies [T] and by Ford and Vilenkin f3l that appeared some years later, ra- 
diation reaction force on moving boundaries attracted the attention of many physicists. 
Due to the movement of the boundary, this topic is also referred to as the dynamical 
Casimir effect (DCE), a name coined by J. Schwinger in his attempt to explain sono- 
luminescence in the early 90s [4J. For a review on this subject see the book by K. A. 
Milton Q and on DCE see Refs. JUlJ!]. 

Though the Casimir force on a unique static plate in vacuum is zero, the fluctu- 
ations of this force are non-vanishing |i9J. Hence, if this plate starts moving with a 
non-zero general acceleration, we expect that a dissipative force proportional to these 
fluctuations appears ifTOl [TTl [T2l . and arguments based on energy conservation lead 
directly to real particle creation. 

Though the DCE already occurs for a unique moving boundary, oscillating cavities 
in parametric resonance with a particular field mode of the corresponding static cavity 
may enhance significantly the particle creation rate lfT3l[T4l[T5l . This effect was studied 
by several authors considering the case of the 1 + 1 ideal cavity ifTSlfTSl . The 3 + 1 
case was also investigated, and different geometries were taken into account, among 
them parallel plane plates |[l3l[ll[16l, cyUndrical [ITJ, and spherical 1, 1 8l [191 |20l IITI 
cavities. The nonideal case was also considered in Refs. Il22l|23l . 

Concerning the static scenario, T. H. Boyer (24] was the first to consider the case of 
mixed boundary conditions (BCs). He demonstrated that the electromagnetic Casimir 
force between a perfectly conducting plate and an infinitely permeable one is repulsive 
rather than attractive. An analogous result was also obtained in the case of a scalar field 
confined within two parallel plates [25. ,26. ,27 J and submitted to a Dirichlet BC at one 
plate and to a Neumann BC at the other. 

The measurement of a repulsive Casimir effect has been persued for many years 
and has finally been achieved very recently by Munday, Capasso and Parsegian |.31i 
in a remarkable experiment involving three distinct media, with appropriate values for 
their permitivity. Although the set up used by these authors in their experiment on 
repulsive Casimir effect is quite different from the two-plate set up made of a perfectly 
conducting plate and an infinitely permeable one, we may learn many things studying 
the DCE with such mixed BCs. Further, though mixed BCs are relatively common in 
the study of the static Casimir effect 1.25. ,26. ,27] l28l |29l [30l . and also in correlated 
topics of Cavity QED 1321 [33l [Ml ESI, the same does not occur for the DCE. In fact, 
as far as we know, the DCE in a 1 + 1 dimensional resonant cavity with mixed BCs 
was considered only very recently, in Ref(s). [36. .37.1 . However, mixed BCs have 
never been considered in the study of DCE for different geometries as, for instance, in 
concentric (and oscillating) spherical shells. 

In a recent paper pTl, the DCE was examined for a massless scalar field submitted 
to Dirichlet BCs at two concentric spherical shells, each of them posessing a time- 
dependent radius. A general expression for the average number of created particles was 
derived for arbitrary laws of radial motions of the spherical shells. Such an expression 
was thus applied to breathing modes of the concentric shells: when only one of the 
shells oscillates and when both shells oscillate in or out of phase. The purpose of 
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this paper is to complement the previous one fTP\ by considering mixed BCs. We 
observe that the field modes associated with mixed BCs differs from that following 
from Dirichlet-Dirichlet BCs. Considering again an oscillatory motion of the shells, 
we identify all the resonances within mixed BCs and derive the expression for the 
associated particle creation rate. Then, performing a numerical analysis we compare 
our results with those presented in Ref. II2TI . For convenience, we shall assume that 
the spherical shell which imposes a Neumann BC to the field is at rest, while the other 
one, which imposes a Dirichlet BC to the field is in arbitrary motion. However, we 
shall consider two situations: in one of them, the inner shell is at rest while the outer 
one is in arbitrary motion and in the other one, the reverse occurs, namely, the outer 
shell is at rest while the inner one is in arbitrary motion. Comparisons of our results 
with those involving only Dirichlet-Dirichlet BCs will be presented graphically. This 
paper is organized as follows: in Section 2 we briefly summarize the main steps of 
the method emplyed to the case where only Dirichlet BCs were considered; in Section 
3 we apply this method to the case of mixed BCs and obtain our general formulas; 
in Section 4, with the purpose of obtaining explicit results for the average number of 
created particles, we choose a particular motion for the oscillating shell and Section 5 
is left for the concluding remarks. 



2 Dirichlet-Dirichlet BCs 

In Ref II2TI the DCE for a massless scalar field confined between two cocentric moving 
shells was considered. This quantum scalar field obeys the Klein-Gordon equation 
□0(r; t) = 0. Besides, this field and its canonical momentum 7r(r; t) = 0(r; t) satisfy 
the equal time commutation relations 

[0(r;i),^(r';t)] = tS{r-r'), 



[0(r;t),0(r';t)] = Kr; t), 7r(r'; t)] - 0. 
The spherical symmetry of the problem leads us to the following solution 
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Fis{r;t) [aims{t) Yim{0,ip) + h.c] , 



ra— — l s — 1 



Fis (r ; t) [aims (t) Yim [e,<p) - h.c] , (2) 



where {Yl„l{^^ ^p)} are the spherical harmonics and the orthonormal radial functions 
satisfy the following differential equation 



1 d (^^ms{r;t)\ , (ujlit) + 



dr 



dr 



Fis{r;t) = {). 



(3) 



Moreover, the operators aims (t) and aj^^ {t) obey the standard commutation re- 
lations 



air,is{t),al„^,^,{t) 



— ^ll'Smm'^s 
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[a/ms(i), arm's' W] = "Ls 4m's' W 



= 0. 



(4) 



Through the time derivative of Eqs. (|2|i, together with the Klein-Gordon equation 
and the canonical momentum formula, we obtain the time evolution for the operators 



aims{t) ~ —iuJls{t)aims{t) + fJ-l\ss']{t)aims' jt) 

s' 



(5) 



where the functions /i;(ss')(i) = [f^iss'it) + ^iis'sit)]/2 and ^ll[ss']{t) = [^J-lss'it) 
— iiis'sit)]/^ are the symmetric and antisymmetric parts, respectively, of the time- 
dependent coefficient 



2uJu[t) 
+ {l-6ss')\l 



rait) 



it) 



r^Fis'{r;t)Fisir;t)dr. 



(6) 



As demonstrated in Ref. |21 1, by comparing Eq. (|5]l with the Heisenberg equation 
of motion dims{t) — i [Hef f{t)iO-ims{'t)] and assuming the most general quadratic 
form of an effective Hamiltonian, we derive 



i 
2 
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a 



l{ — ra)s'j Ims 



l,m,s,s 



(7) 



The evolution of the density operator is computed through the relation p{t) = i [Hef / (t) , p{t)], 
with the aid of an iterative procedure up to second order approximation in the velocity 
of the cavity boundaries, i.e., fi{t), ro{t) <C c. The derivation of the average number 
of particles created in a particular mode — labeled by the quantum numbers (I, m, s) 



— is thus given by Uims{t) = Tr p{t)a]^^{{))aims{Q) 
state p(0) = I {0}) I {0}) it follows that 



, and for an initial vacuum 



Mms(i)=X! / dii^/(s/s)(ti) exp{i + rj;s(ti)]} 



(8) 



with Vtisit) = dti ujisih). 

The number of created particles Nims{t) depends on the radial function Fis{r; t) 
through /ii(s's) (t). The solution of Eq. ^ is given by a linear combination of spherical 
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Bessel functions of the first ( ji) and second (rt;) kind, such that the Dirichlet BC appUed 
to the inner shell leads to the relation 

Fis{r;t) = Ar,,[jj^,(t)r)n,^,(t)r,(t)) - j,^,(t)r,(t))n,^,(t)r)], (9) 

whereas that on the outer shell results in the transcendental equation 

i/^.Wro(t))7i;^,(<)r,(<)) - j,^,(t)r,(t))7i,(c^,,(t)ro(t)) = 0. (10) 

In Fig. 1 we present a map of the solutions of Eq. (flOl i for some values of the numbers 
/ and s. As it was noted in ||2T1 . for the case I = 0, the frequencies are equally spaced. 
This fact does not occur for the case I ^ 0. However, when both radii of the shells 
are much larger than the separation between them, i.e., Vi it) » Tq (t) — (t), the 
solutions for all values of I approach the solution for the onedimensional case, so that 

UJls sit/ {To (t) - Ti (t)). 



o 
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Figure 1: Map of the solutions of the transcendental equation ( fTOl l. The colors corre- 
spond to different values of the number I. The solid, dashed, and dotted lines corre- 
spond to s = 1, s = 2, and s = 3, respectively. 



3 Mixed boundary conditions 

It is important to emphasize that the expression for the average number of created 
particles, derived in Eq does not depend on the character of the BCs. Thus it can be 
applied even for mixed BCs as in the present work, where we assume that the massaless 
scalar field satisfies the Neumann BC at a fixed spherical shell and a Dmchlet BC at a 
second concentric spherical shell whose radius has an arbitrary time dependence, 

a.0(r,t)|,=,„ =0 and (r, = 0, (11) 

where the index a {(3) is related to the static (moving) shell. However, different expres- 
sions come up for Fis (r; t) and ujis (t), as compared to those in Ref. ( fTTI i. 

As already noted in the previous section, the general solutions to Eq. ^ are lin- 
ear combinations of spherical Bessel functions, but the mixed BC leads to a different 
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solution. The assumption of Neumann BC for the field at the static shell leads to the 
following expression for the radial functions 



f d d \ 

Fis{r;t) = Nis ji {uJis{t)r) - — ni {ujis{t)ra) - m {ujis{t)r) - — ji {uJis{t)ra) , 

(12) 

and the subsequent assumption of Dirichlet BC on the field at the moving shell leads 
to a frequency discretization 

d d 

- — ji{uJis{t)ra) ni{ujis{t)rf3{t)) = ji(uJis{t)rp{t)) - — ni{ujis{t)ra) ■ (13) 

In Figs. 2 and 3 we show the maps of the numerical solutions of the transcendental 
equation Eq. (T3[ for some values of I and s. As we can see, the map of w/s (t) is very 
sensitive to the BCs. For Dirichlet-Dirichlet BCs (Fig. 1) the frequencies LjJos{t) are 
equally spaced, a situation that does not occur for mixed BCs. We also note that the 
map of ujis (t) turns out to be entirely different when considering the Dirichlet BC in 
the outer shell and Neumann BC in the inner one (Fig. 2) or oppositely, with Neumann 
BC in the outer shell and Dirichlet BC in the inner one (Fig. 3). 
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Figure 2: Map of the solutions of the transcendental equation ( fT3T l with < r^j. The 
colors correspond to different values of the number I. The solid, dashed, and dotted 
lines correspond to s = 1, s = 2, and s = 3, respectively.. 



However, there are also some similarities between the results derived from mixed 
BCs and those for Dirichlet-Dirichlet BCs, in Fig. 1; it can be directly verified that the 
solutions for ujis{t), coming from Eq. (13), approach that for the onedimensional case 
(ojis — *■ (s — 1/2) it/ Itq, — {t)\) when both radii of the shells are much larger than 
the separation between them. 

A comment is in order here: we note that the BCs must be imposed in the instan- 
taneously co-moving Lorentz frame, where the boundaries are momentarily at rest. If 
the Neumann BC was to be imposed on the moving boundary, we should have used 
the appropriate Lorentz transformation to write the fields in the inertial frame of the 
laboratory as follows 

ar'</'(r',i) I. =^ R + f^(i)aa0(r,t)| =0. (14) 
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Figure 3: Map of the solutions of the transcendental equation ( fT3] l with > r^. The 
colors correspond to different values of the number I. The solid, dashed, and dotted 
lines correspond to s = 1, s = 2, and s = 3, respectively. 

In that case, the time derivative in Eq. ( fT4b invalidates the expansion used in Eq. (|2]l, 
and as a consequence, also in Eq. This fact demands a different formal develop- 
ment for the computation of the required particle creation rate. For that reason, in the 
present work we treat only the case where the Neumann BC is imposed on a spherical 
shell at rest, leaving aside the breathing modes analyzed in Ref. ITTll . when both shells 
oscillate in or out of phase. 



4 Numerical estimatives 

In this section, in order to obtain explicit results, we will consider an specific motion 
for the spherical shell that imposes on the field the Dirichlet BC. A typical situation 
consists of an oscillation that starts at some instant, has a sinusoidal behaviour with an 
angular frequency w and a small amplitude and then stops at some later instant. We, 
then, assume that the radius of the moving shell has the following law of motion 

r(i{t) = r(i(\ + (.sni{wt)), (15) 

with e -C 1. In the following we also assume that the cavity mirror oscillates only 
during a finite time interval T, then stopping suddenly its motion. 

Substituting Eqs. (|6]l and (fTsT i into Eq. dHJ and making a power series expansion 
with respect to the small parameter e, we obtain 

Mi^s^(^^ Y^\Ciiss')!isA^\T)^ , (16) 

where, after defining ujiss' = ijJisi^) + ^^/s'(0), the coefficient C/ss' and function 
jiss' (07; T) are given by 

&^^\i[w - ijJiss')T\-\ exp[-j(tx7 + u;;5s/)r] - 1 
Jlss' [T^] -L ) ^ — ■ — , (17) 

1{W - LUlss') -L I (W + UJlss' ) -L 

and 

1 diuisiO) 



1 



- rp{l-5ss')\ 77^:/ dr r F/^ r; . (18) 

Note that fiss' {vj; T) is an oscillating function of T, except when the mirror os- 
cillating frequency satisfies the resonance condition, namely, m = loiss' ■ In this case 
fiss' i^^iss' ',T) — 1, and the number of created particles turns to be the following 
quadratic function of T, 

lim Aftms^ f!^^!!f£l\ \\ (19) 

This result is valid only under the short-time approximation eujiss'T << 1, since we 
have disregard terms proportional to (euJiss'T)^ with n > 3. Moreover, Eqs. ( fTSI l to 
( fT9] l are valid either for Dirichlet-Dirichlet BCs or mixed BCs since they were derived 
using the equation of motion ( flST l and equations dTJ - dHJ which are independent from 
the BCs. 

To study the behavior of our result in Eq. il9[ . we plot in Fig. 4 the expression 
J^ims/ (eroT)^ as a function of njTT/{ro — ri) for some values of I and s, under reso- 
nance conditions zu = ujiss', setting ro = 2ri. Both letters in the legend indicate the 
BCs on the inner and the outer shells, respectively: for example, D (D) means Dirich- 
let BC on a static (moving) shell, whereas N means Neumann BC on a static shell. As 
we can see, both the intensity and position of the resonances change in a non-trivial 
way with the BC. The case of a moving outer shell with the field satisfying Dirichlet- 
Dirichlet BCs exhibits higher resonance intensities, while the case of a moving inner 
shell with the field submitted to mixed BCs leads to lower resonance frequencys. 

In the limit Tq — r^, we can use the Bessel asymptotic forms for large argu- 
ments to derive an analytical expression for the average number of created particles in 
a particular mode. For the case where the field is submitted to Dirichlet-Dirichlet BCs, 
we obtain 

L^iss' T ^ (20) 

and 

lim Mhns — > -. rrs's- (21) 

Analogously, for the field submitted to mixed BCs ( fTTI ). we have 

(g + g'-l)7r 

^Iss' — > i i (^^) 

\ra - rp] 



and 



lim Ni„,s ^ ''""'7' , ,4 {2s' - l)(2g - 1). (23) 



Expressions ( |20] | and (ISTT l correspond to the results for the 1 + 1 DCE under Dirichlet- 
Dirichlet BCs derived in Refs. lfT3l[T5l l6ll. whereas Eqs. (|22] | and (|23] | correspond to 
the results under mixed BCs presented in Il25ll26ll27l . These similarities can be related 
to the fact that the limit r; ^ Tq — is akin to the plane geometry. 



8 



I 0.2- 



0.1- 







• DD 


■ DN 


♦ DD 


* ND 




♦ < 




A 



T r 



^ 0.3n 



0.2- 



0.1- 



Eg 0.3n 



0.2- 



"T" 



— r 



' * * 

~i 1 1 



10 



WTT/{ro ~ ri) 



Figure 4: Plot of Afims / (e'coT)'^ as a function of -cutt/ {to — Vi), in the resonance con- 
dition for a few values of I and s. We have considered Dirichlet-Dirichlet and Mixed 
BCs. On the legend (top-right of the figure), the letter on the left indicates the BC 
imposed on the field at the inner shell and the letter to the right, indicates the BC im- 
posed on the field at the outer shell. D means Dirichlet BC and static shell, D means 
Dirichlet BC and moving shell, whereas N means Neumann BC and static shell. We 
have set Tq = 2ri. 
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5 Concluding remarks 



In this paper we have investigated the dynamical Casimir effect for a massless scalar 
field within two concentric spherical shells considereing mixed boundary conditions. 
We have thus complemented some previous results presented in Ref. [21 J where the 
massless scalar field was assumed to satisfy only Dirichlet BC in both shells. We have 
analyzed the real particle creation phenomenon for the case where only one of the 
shells is allowed to move with an arbitrary law of motion for its radius. In addition, 
the Dirichlet BC was imposed on the moving shell while the Neumann's was assumed 
on the static one. However, in our discussion, the moving shell could be the inner 
shell or the outer one as well. In order to get some numerical estimatives, and with 
the purpose of comparing our results with those obtained in Ref. ( l|2ll ). we chose a 
particular, but very typical, oscillating motion for the moving shell, in which it starts 
moving at a certain instant, oscillates with a given frequency and then stops suddenly 
its motion. Considering this particular situation, we have identified the resonance con- 
ditions where the number of created particles is more appreciable. A direct inspection 
in our graphs (see Fig. (4)), allows us to make some conclusions: for both cases of 
Dirichlet-Dirichlet BCs or mixed BCs, we see that every time the moving shell is the 
outer one the average number of created particles is greater than the corresponding 
cases where the inner shell is in motion (by a factor of the order of 4) . This can 
be unsderstood simply recalling that the dynamical Casimir effect increases with the 
area of the moving surface. In other words, the dissipative force that acts on the mov- 
ing boundary, responsible for converting mechanical energy into field energy (real field 
quanta) increases with the area with the moving boundary. Another interesting result 
that can be extracted from our calculations is the fact that the case with mixed BCs 
presents lower resonance frequencies than that with Dirichlet-Dirichlet BCs. This fea- 
ture can be useful for further experimental investigations of particle creation within the 
context of the dynamical Casimir effect, since it makes easier to access the parametric 
amplification regime of paiticle creation. 
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